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Abstract 

Berezinskii-Kosterlitz-Thouless (BKT) transition, the transition of 
the 2D sine-Gordon model, plays an important role in the low dimen- 
sional physics. We relate the operator content of the BKT transition to 
that of the SU(2) Wess-Zumino-Witten model, using twisted bound- 
ary conditions. With this method, in order to determine the BKT 
critical point, we can use the level crossing of the lower excitations 
than the periodic boundary case, thus the convergence to the tran- 
sition point is highly improved. Then we verify the efficiency of this 
method by applying to the S=l,2 spin chains. 

PACS numbers: 75.40. Cx, 05.70.Fh, ll.10.Hi, 75.10.Jm 
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1 Introduction 



Berezinskii-Kosterlitz-Thouless (BKT) transition ^ |3| plays an important 
role in 2D classical and one- dimensional (ID) quantum systems, such as the 
2D XY spin model, 2D helium film, 2D superconducting film, roughening 
transition, ID quantum spin models and ID electron systems. Furthermore, 
BKT transition is a typical instability of the Tomonaga-Luttinger liquid, such 
as the Mott transition in ID metal [[|. 

It was a difficult problem to determine numerically the critical point and 
the universality class of the BKT transition, because of logarithmic correc- 
tions and the slow divergence of the inverse correlation length (or the energy 
gap) ||. Moreover, it is hard to distinguish the BKT transition from the 
second order transition, by using the conventional finite-size-scaling method 
H [7]. In fact, for the exactly solvable S=l/2 XXZ spin chain, where the 
transition at A = 1 is of the BKT type, the finite-size scaling method leads 
to false conclusions ||. This problem was successfully resolved by the level 
spectroscopy |J, based on the renormalization group calculation and the 
SU(2)/Z2 symmetry in the BKT transition. 

The SU(2) symmetry inherent in the BKT transition was pointed out 
by Halpern |TIJ and Banks et al. [|TTJ. They showed the equivalence between 
the SU(2) massless Thirring model and the theory of the bosons consisting 
of a free field plus a f3 2 = 8tc sine-Gordon model, which corresponds to the 
BKT line. Another approach was proposed by Ginsparg [0, based on the 



c = 1 conformal field theory: modding out the SU(2) symmetry model by Z 2 
symmetry gives the structure of the BKT multicritical point. 

In this paper we directly relate the SU{2)/Z 2 structure of the oper- 
ator content of the BKT transition to that of the k = 1 SU(2) Wess- 
Zumino-Witten (WZW) model, using twisted boundary conditions. With 
this method, in order to determine the BKT critical point, we can use the 
level crossing of the lower excitations than the periodic boundary case, thus 
the convergence to the transition point is highly improved. We then apply 
the method to the S=l,2 spin chains. 

2 BKT transition and twisted boundary con- 
dition 
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2.1 On the Gaussian line 



We consider the 2D Gaussian model defined as the Lagrangian 

We compactify on a circle as = + 2n/ y/2. We introduce the dual field 
9 to defined as 

9,0 = -d v {iK0), d y <j> = d x (iK9), (2) 

which has the periodic nature 9 = 9 + 2ix / y/2. 
The vertex operator defined as 

O m ,n ='■ exp(imy/2(f)) exp(iny/26) :, (3) 

has a scaling dimension x mj „ and a conformal spin s mjTl 

x m ,n( K ) = \ {rn 2 K + ^ , S m>n = ffffl. (4) 
There are the other type fields, i.e., the current fields 

90, (5) 



y/K y/K 
which have x = 1, s = ±1, and the marginal field 

jUf = -1^090 = -l ((9,0) 2 + (dy<P) 2 ) (6) 

which has x = 2, s = 0. There exist also descendant fields of them. 

About the symmetry, the Gaussian model ([[]) is invariant under — ► 
4>+const., 9 — > 9+const., which means the £7(1) x 17(1) continuous symmetry. 
There are also the discrete Z 2 symmetries, T : (z, 0, 9) — > (z, —0, —0) and 
C : (z, 0, #) — ► (z, 0, — 9) (we can also define the redundant symmetry P = 
CT : (z, 0, 9) — > (z, —(f), 9)). The Gaussian model is invariant under the dual 
transformation K <-> 1/K, <f> «-> 0. The self dual point if = 1 is nothing 
but the k — 1 577(2) x SU(2) Wess-Zumino-Witten model, whose symmetry 
structure is apparent in the conformal dimensions of the vertex, the current, 
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the marginal operators and their descendants. At the point K = 1, the 
number of the marginal operator (x = 2, s = 0) is 9. 

The K = 4 is the BKT transition point, where the number of the marginal 
operator (x = 2, s = 0) is 5. The correspondence of the scaling dimensions 
for the WZW point (K = 1) with those for the BKT point (K = 4) is given 
by 

£2m,n(l) = ^m,2n(4), (7) 

and 

£m,2n(l) = £n,2m(4). (8) 

In short, there is a partial correspondence of the operator content of these 
two models, reflecting the Z 2 symmetry difference. 

Fortunately, introducing a half-integer magnetic charge, one can obtain 
the full correspondence of the K = 1 and K = 4 Gaussian models. With 
twisted boundary conditions on a cylinder, or the cut from to infinity on 
the complex plane, one can set the arbitrary charges at ±oo (on the cylinder) 
or at 0, oo (on the complex plane) |13|, [14], [L5| (see Appendix A). The relation 
between the twisted boundary condition and the change of the charges has 



been observed in the Bethe-Ansatz solvable S=l/2 XXZ spin chain [16 
Therefore, with the twist angle $ = tt, one can introduce a half-magnetic 
charge. In that case, we can set 



x, 



T m B n c (K)=x m+1/2 , n (K). (9) 



Then, we obtain 



x 2m+ltn (l) = x 1 m f n (A), (10) 

and 

X m,2n(l) = s n,2m+l(4), X m2n+ i(^) = ^n,2m+l(4) ; (H) 

so that we can see the correspondence of all the operator content of K — 1 
and K = 4 points, considering twisted boundary conditions. 

2.2 Renormalization 

Next we introduce the sine-Gordon type interaction, that is, 
Ci = : cos V84> : (for K = 1), : cos v 7 ^ : (for K — 4), (12) 



4 



(where a is a short- distance (ultraviolet) cut-off). Here we assume that there 
is no cos y/2(J) interaction term for the K — 1 case, by the symmetry reason 
(discrete symmetry — > + 7r/\/2 which corresponds to "translation by one 
site" in the half-odd integer spin chain [pjj) or adjusting parameters. With 
the simple transformation —>■ 20, 9 — > 9/2 or m — > 2m, n — > n/2, the sine- 
Gordon model for K = 4 becomes equivalent to the sine-Gordon model for 
K = 1. 

About the symmetry, the U(l) symmetry of is explicitly broken to the 
discrete symmetry 0^0 + 2tt / 'y/8(K = 1) or — > + 2tt/V2(K = 4). 
Furthermore, the Lagrangian ([D) is invariant under — > + 7r/v8, j/^ — > 
-^(■fif = 1) or -> + vr/v/2, ^ -> -j/^tf = 4). 

Under a change of cut-off a — ► e l a, the renormalization group equations 
for the sine-Gordon model are || 

dy (l) 



dl 
dl 



-y%i) 

-V&)Vo{l) (13) 



where K = 1 + y /2(around K = 1) or if = 4(1 + y /2) (K = 4). For 
the finite system, / is related to the system size L by I = log(L/a). There 
are three critical lines; y^ — corresponding to the Gaussian fixed line, and 
2/0 = ±?/o (yo > 0) corresponding to the BKT lines. On the BKT lines, the 
couplings behave as yo(l) = ±2/^(0 = + l/i/o(0)) = l/log(L/L ). In the 
region between the two BKT lines, all the points will be renormalized to the 
Gaussian fixed line, so they are massless. The other region is massive, except 
on the Gaussian fixed line. 

Although the renormalization flow is the same between the K = 1 and 
the K = 4 cases, the operator content is different, and in general the corre- 
spondences (H), ( |ITD are not satisfied, since the duality relation K <-> 1/K 
does not hold on the BKT lines y$ = ±j/ . Fortunately, the correspondences 
(0), ( |T0D remain correct after the renormalization, since the operator product 
expansion (OPE) structure of both cases is the same, for example, 

K = 1 K = A 

(cosv / 20cosV / 80cosV / 20) <-» (cos0/-\/2cos V / 20cos0/v / 2) 

(sin v^20 cos a/80 sin V2(j)) <-> (sin^/v^cosv^sin^/v^) 

(cos cos Vo0) <-> (cos V^M cos \/20) (14) 
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which can be seen by the mapping — > 20, 9^9/2 (similarly, higher OPE's 
are the same) . The renormalized scaling dimensions are determined by these 
OPE |18|, |9[ (see Appendix B). Therefore, using twisted boundary conditions, 
we can explicitly see the SU(2) structure inherent in the BKT transition. 

In Table 1 we summarize the operator content of the WZW (K = 1) 
model and that of the BKT (K = 4) model, with renormalized scaling di- 
mensions. Note that at the BKT critical line (?/</,(/) = 2/o(0)> there are de- 
generacies of the excitations corresponding to x pl,p2 ,x p3 or x m0 ,x m3 ,x m7,m8 
or x m2 ,x m5,m6 , respectively, reflecting the SU(2) (SU(2)/Z2) symmetry (for 
the y^l) = —yo{l) branch, the role of operators x p0 «-> x p3 , x m5,m6 <-> x m7,m8 
interchanges). Thus the level crossing of them can be used to determine the 
BKT critical point. 

In practical systems, there are corrections from the descendant fields of 
the identity operator 1. The most important irrelevant fields in them are 
L_ 2 ^-2l, (L 2 _ 2 + L 2 _ 2 )l with scaling dimension x = 4 |0j. With the 
twisted boundary conditions, we can use the level crossing of the lower ex- 
citations than that with only the periodic boundary conditions. Since the 
amplitude of the corrections from the irrelevant field becomes smaller with 



the lower excitations [ 20 1 , thus the convergence to the transition point is 
highly improved. 



3 Physical examples 
3.1 S=l XXZ chain 

As a physical example, let us consider the S=l XXZ spin chain, described 
by the Hamiltonian: 

L 1 

H = hjj+i, h jtj+1 = ~(SfS j+1 + Sj Sf +1 ) + ASjSj +1 , (15) 

where we assume L is even and periodic boundary conditions. This model 
is invariant under spin rotation around the z-axis, translation(Tft : Sj ,y,z — > 
Sj+{ z ), space inversion (P : Sj ,y,z — > S^ ,y f +1 ), spin reversal (T : Sj — > 
Sj,Sf -> SJ). Therefore, energy eigenstates are characterized by the z- 
component of the total spin (S? = J2S Z ), wavenumber (q = 2nk/L), parity 
(P = ±1), spin reversal (T = ±1). Note that with the unitary transformation 
exp(7rz jSj )) spin operators change as Sj — > Sj,Sf — > (—iysf, thus the 
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momentum of the excitation in the odd changes q — > q + it, and the sign 
of the XY term in flT5|) reverses. 



About the spin S XXZ chain, Haldane [21| discussed that for the integer 
S 1 case, there are two transition points A c i < 1 < A C 2, where the XY-Haldane 
transition A c i is of the BKT type, and the Haldane-Neel transition point A C 2 
is of the 2D Ising universality class, in contrast to the half-odd integer S case 
where only one transition point exists at A = 1 of the SU(2) WZW type. 

Numerically Botet and Jullien estimated A c i ^0.1 [p2|, Sakai and Taka- 



hashi A c i = -0.01 ±0.03 [|3|, and Yajima and Takahashi A = 0.069 ±0.003 
|24|| . Recently, the authors found A c i = and checked the universality class 
|25| using the level spectroscopy |J. We found the exact degeneracy (at least 
within numerical accuracy) between the = ±4, q = 0, P = 1 and the 
Sf. = 0, q = 0, P = T = 1 excitations. 

These results are obtained under periodic boundary conditions. To con- 
sider a twisted boundary condition, in ( |15| ) we replace the boundary term 
between the L-th site and the 1st site, 

h L ,i = \(S+S^ expH$) + 5 L 5 1 + exp(z$)) + ASffi. (16) 



However, with this boundary condition the Hamiltonian (|T5| ) is not transla- 
tional invariant. The translational invariance can be restored, still maintain- 
ing the total angle $, by twisting all neighboring bonds in the chain by an 
angle $/L with the next unitary transformation 

/ $ 1 

u* = ex P fi-x;(j--)s; 

UzSfU* 1 = Sfexp(±i(j -1/2)$/ L), U<s,S*U$ l = S*, (17) 

(in [/$ we use (j — l/2)Sj for the compatibility of the definition of P). Thus, 
we obtain 

H t(®) = E Mj+i. h h+i = \{SpJ +l e-^l L + SjS+ +l e^/ L ) + AS*S* +1 . 

(18) 

Under this boundary condition, we had better modify the definition of the 
translational operator T l R = exp(i$S^/L)TR (see Appendix C). Although 
there are no discrete symmetries like P, T in general for the special 
angle $ = n, the system (|i~8l) is invariant under the discrete symmetries 
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U 27r P, U 27T T. Moreover, it was shown that U 27r P, U 2n T are the good quantum 
numbers |26| characterizing the generalized Z 2 x Z 2 symmetries [27[] (see also 



Appendix C). Thus the twist angle $ = ir has the special meaning in addi- 
tion to the one discussed in section 2. Hereafter we call the twisted boundary 
condition with $ = 7r as TBC and the periodic boundary conditions as PBC. 
The energy eigenvalues E n are related to the scaling dimension x n as 

E n (L)-E g (L) = 2 -^, (19) 



where v is the spin wave velocity |28[ . And the conformal anomaly c is related 
to the ground state energy O, p9| 



TT11C 

E g {L) = e 9 L - — . (20) 

In Fig.l we show the excitations corresponding to the scaling dimension 
x = 1/2. The excitation Sj> — 0, q — 0,U 2w P = U 27T T = —1 under TBC is 
exactly degenerate with excitations Sj> = ±2, q — 0, P — 1 under PBC at 
A = 0. Fig.2 shows the excitations corresponding to the scaling dimension 
x = 2. At A = 0, the excitations S% = ±2, q = 0, U 2w P = -1 under TBC 
are exactly degenerate with excitations Sj. = 0, q = 0, P — T = 1 and 

= ±4, q = 0, P = 1 under PBC, whereas the excitations = ±2, q = 
0, U 27T P = 1 under TBC are exactly degenerate with the excitation = 
0,q = 0,P = T = — 1 under PBC. Fig. 3 shows the excitations with scaling 
dimension x = 1, = 0, q = ±2ir/L with PBC, S Z T = ±2,q= ±2n/L with 
TBC. They are exactly degenerate at A = 0. 

Next we discuss the universality relations. The conformal anomaly is 
estimated c = 1.000 in p^fl . We can eliminate the logarithmic corrections 
by taking the appropriate average, for example, (x p0 + 2x pl + x p3 ) as can 
be read from Table 1. We show the size dependence of (x p0 + 2x pl + x p3 )/A 
in Fig. 4; the remaining size dependence is mainly explained by the x = 4 
irrelevant field. From Table 1, in the neighborhood of the BKT transition, 
there appear terms linear in the distance t from the BKT line. The ratios 
of x p3 - x pl , x ml - x m3 , x m5 - x m2 } x m0 - x m3 , are -1/2 : -1 : 1 : -4/3 
close to the BKT line. In Table 2, we show the coefficients linear in t of the 
eigenvalues corresponding with x p3 — x pl , x ml — x m3 , x m5 — x m2 , x m0 — x m3 . 
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3.2 S=2 XXZ chain 



Next we consider the S=2 XXZ spin chain, described by the Hamiltonian: 

H = J2(S*S* +1 + S]S] +1 + ASjSj +1 + D{S*) 2 ). (21) 

For the isotropic case (A = 1, D — 0), several studies have been done in 
the relation to the Haldane's conjecture[|2l| |30|, [H], |32], ^3], [34], |35], 



Although the estimated values of the Haldane gap are widely ranged, they 
are common on the existence and the smallness of it. For the smallness of the 
Haldane gap, we can expect that the transition points between the Haldane 
gap and the XY phases 1 — A c (D = 0) and D c \ (A = 1) are very small, 
compared with the S = 1 case[]3"8"|, Here we only consider on the line 
D — 0, and on the line A = (D > 0), to estimate the XY-Haldane and the 
XY-large D transition points and to determine the universality class. 

First we consider on the D = line. Figure |5] shows the excitation energies 
with S% = ±2, P = 1 under PBC, and with Sf, = under TBC for L — 12 
systems. We can see a level crossing which corresponds to the transition 
point between the XY and the S = 2 Haldane phases, as is expected from 
Table 1. The size dependence of this crossing point is shown in Fig. |6|. 
The extrapolated value to the thermodynamic limit is A c = 0.966. The 
conformal anomaly number of this point is estimated as c = 1.16. To check 
the universality class, we show the size dependence of the combination of the 
scaling dimension as (x p0 + 3x pl )/4 in Fig. |7]. We extrapolate the L — > oo 
as A C (L) = A c + ai/L 2 for L = 10, 12 and A C (L) = A c + ajL 2 + a 2 /L 4 for 
L — 8, 10, 12, and the obtained values are 0.487 and 0.491 respectively. 

Next we show the result on the A = 1 (D > 0) line. Figure [| shows 
the excitation energies with = ±2, P = 1 under PBC, and with Sj> — 
under TBC for L = 10 systems. We can recognize the three regions as the 
S = 2 Haldane phase < D < D c i, the massless XY phase D c \ < D < D c2 , 
and the large D phase D c2 < D. The size dependence of crossing points 
is shown in Fig. |9| and 10, and the estimated values are D c \ = 0.043 and 



D c2 = 2.39. The estimated conformal anomaly numbers are c = 1.16 and 
0.998 for D c i and D c2 respectively. The transition point between the S — 2 
Haldane gap and the XY phases (D&) is consistent to the previously obtained 
values by Schollwock and Jolicceur [pi] (D c i = 0.004(2)), but the transition 
point between the XY and the large-D phases (D c2 ) deviates from their value 
(D c2 pa 3). Figure O and O shows the size dependence of (x p0 + 3x pl )/4. 



The extrapolated values for D cl are 0.489 for L = 10, 12 and 0.493 for L 
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8, 10, 12. For D c2 the values are 0.5005 for L = 10, 12 and 0.5006 for L = 
8,10,12. 

The conformal anomaly numbers and the scaling dimensions for the crit- 
ical points (A, D) = (A c , 0) and (1, D c i) somewhat deviate from the ideal 
value c = 1 and x = 1/2. We think that at these points the BKT transition 
points are very close to the 2D Ising transition point between the 5* = 2 Hal- 
dane and the antiferromagnetic phases@], so that the crossover effect is very 
large for these points. In fact, at the BKT transition point (A, D) = (1, D c2 ) 
where the 2D Ising transition point is far, these values are consistent to the 
ideal values. 

In the XY phase on the A = 1 line, we find that the energy levels with 
= 0, U 2n P = U 2w T = 1 and with = 0, U 2n P = U 2n T = —1 cross two 
times (see Fig.[l3|). These two points correspond to the Gaussian fixed points 
[ 15 1, so that there are two Gaussian fixed lines in the whole phase diagram 
of (|2~ip. This may be the "indirect" evidence of the intermediate-/} phase, 



predicted by Oshikawa |27 [. To see the intermediate-/} phase, Oshikawa, 



Yamanaka, and Miyashita|40| studied the line A = 1 with the quantum 
Monte Carlo method. But we cannot find that phase, and we think that the 
intermediate-/} phase needs more large values of A (> 1). 

Lastly we remark the following thing. From the large S mapping onto the 
anisotropic nonlinear sigma model by Haldane||21||, 1 — A + D is proportional 
to the anisotropy of it. Numerical values 1 — A c = 0.034 and D c i = 0.043 
are comparable, thus it is consistent to Haldane's arguments. 



4 Conclusion 

Physical properties of the BKT transition, including the renormalization 
group properties and the SU(2) or the SU{2)/Z 2 symmetry, have been well 
investigated in the field theory. However, the mapping from the various 
models to the field theoretical models, such as the sine-Gordon model or the 
Wess-Zumino-Witten model, is not simple, not quantitatively correct. Fortu- 
nately, the symmetry structure and the sum rule at the BKT phase transition 
point survive after the mapping. Therefore, using these properties, we can 
determine the BKT critical point and the universality class. 

In numerical calculations, SU(2) symmetry has been used to determine 
the BKT-type critical line in |41|| , SU(2)/Z 2 symmetry has been used in 



and the sum rule to eliminate logarithmic corrections in [42]. Recently, 
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introducing the half magnetic charges by twisting boundary condition, one 
of the authors has developed a method to determine the Gaussian fixed 
line and its universality class for the non-integrable models [DJ. In this 
paper, using the twisted boundary condition, we explicitly relate the energy 
eigenvalue structure to the SU(2) symmetry of the BKT transition. At the 
same time, we can improve the convergence of the physical quantities to the 
thermodynamic limit, comparing with the original level spectroscopy ||. 

Our method is applicable not only to the quantum problems, but also to 
the classical models, treating the eigenvalue structure of the transfer matrix. 
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A Free boson on the complex plane 

We shortly review the free boson theory on the complex plane ]|3) and the 
relation between the half-odd magnetic monopole and twisted boundary con- 
ditions [14|. The equation of motion ddcf) = for ([J) allows the chiral de- 
composition 

4>{z,M) = ^{<p{z)+q>(z)), (22) 

where we introduce the two independent complex coordinates z = x + iy,z = 
x — iy and use the notation d = (d x — id y )/2, d = (d x + id y )/2. Then the 
action can be written 

The two point functions are 

(<p(z)<p(w)) = - \og(z - w), ((p{z)(p{w)) = - \og(z - w). (24) 

Chiral current: 
Here we introduce the U(l) chiral current as 

J(z) =idip(z). (25) 

The chiral current has a leading short distance expansion 

J(z)J(w)= 1 +•••, (26) 
[z — w) z 



inferred by taking two derivatives of We introduce the mode expansion 
of the current 

J{z) = ^an, a n =i ^z n J(z). (27) 

Using the short distance expansion (^6]) and the radial quantization, we ob- 
tain the commutation relation 

[a m ,a n ] = m8m+nfl, (28) 

which means the U(l) current algebra. Note that «o is the conserved charge 
of the U(l) current. 
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Stress energy tensor: 
From the Noether theorem, the stress-energy tensor is written as 

T(z) = I : J(z)J(z) : . (29) 

Using the short distance expansion (|26|) and Wick's theorem, we obtain OPE 
of the T(z),T(w) 

r( , m „ ) = ^ + «M. + «M +r ^, ,30) 

(z — w) A (z — w) z z — w 

which means that this satisfies the c = 1 Virasoro algebra. Similarly, OPE 
of the T(z),J(w) 

r^t X W \ J(W) dJ(w) , . 

T(z)J(w) = - 1 ; + — ^ + reg. 31 
\z — w) 2 z — w 

means that J(z) is the primary field with the conformal dimension (1,0). The 
mode expansion of T(z) is 

T{z)=Y.z- n - 2 L ni (32) 

n 

where L n can be written 

^ n = o ' a m+n<^-m : ; (33) 
* m 

especially 

oo 

-^0 = + : a ~n®n ■ ■ (34) 

^ 71=1 

Vertex operator: 

Vertex operator is defined as : exp(icap(z)) :. OPE of the vertex operator 
and the current is 

J(z) : e iav{w) := : e ia ^ w) : +re#.. (35) 

OPE of the vertex operator and the stress-energy tensor is 



<T 2 /2 . . r) ■ P i°«p{w) ■ 

T{z) : e %atp ^ := : e* a</,(u,) : + _ 



(2; — w) 2 z — w 
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+ re^., (36) 



therefore the vertex operator is a primary field with conformal dimension 
h = a 2 12. 

Highest weight states: 
When we define the state 



la) = lim : e ia ^ w) : |0), (37) 

w—*0 

it is the highest weight state of the U(l) current algebra, 

ao|a) = a|a), a„|a) = (n > 0), (38) 
from eq. fl3~5|). This is also the highest weight state of the Virasoro algebra 

a 2 

L \a) = —\a),L n \a) = 0(n>0). (39) 

Compactification of the internal space: 
The compactification of the internal space = + 27r/v / 2 restricts the 
eigenvalues a, a of the U(l) x U(l) charges a ,a . First we note that 

<f(z) = q - ia \ogz + % V -z~ n a n , 

<p{z) = q-ia \ogz + i V -z~ n a n , (40) 

where q, q are zero-modes and satisfy the commutation relation 

[q, ao] = [q, a ) = i. (41) 

Rotating around the origin 0, there appears a phase factor ±2m in the func- 
tion log(z)(log(z)). Thus, considering the periodic nature of 0, we obtain 

[J 

a — a = J — n, (n : integer), (42) 

for the requirement of the uniqueness of the vertex operator under the change 
= + 2n/V2, 

a + a = y/2Km, (m : integer). (43) 
In one word, eigenvalues for ao, ao are 



a 



K 1 K 1 

— m + W n, a = \ — m — \ n. 44 

2 \2K 1 V 2 V 2iT V ; 
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Since L ± L are the generators of dilatation and rotation, we identify (a 2 ± 
ct 2 )/2 as the scaling dimension and the conformal spin of the state \a,a). 
Finally, by defining the dual field # to as 

9=^=(<p(z)-<fi(z)), (45) 

9 has a periodicity 2tc/\/2. 

Twisted boundary conditions and Background charge: 
To any conformal operator /(a ,ao), we can associate a twisted operator 

f @ e(a ,a ) = exp(-i(Oq + Qq))f(a , a ) exp(i(6g + Qq)) 

= f(a + e,a + e) (46) 

Here we set 

[k $ 

e = e = fe (47) 



On the one hand, from eqs. (f40|) this means 

^^-^I^^'-^tA^' (48) 

that is, there is a cut for 9. Mapping from the plane to the cylinder by 
w = u + iv = (L/2ir) log z, we obtain 

K $ 1 $ 

^^_ ! __ M ,^ ( , + __„, (49 ) 

that is, twisted boundary conditions for 9 in the v direction. 

On the other hand, the change ([46D , ( flTf ) with eqs. ([42]) , (f43l) means the 
change of the magnetic charges 

m — ► m H , n — > n. (50) 

27T 

This can be interpreted to set a magnetic monopole $/27r at the origin z = 0, 
and a magnetic monopole — <3>/27r at the infinity z = oo. 

Discrete Symmetries 
For the Gaussian model (p3|), besides continuous t/(l) x symmetries, 
there are discrete symmetries 

T : (<p, (p) -> (-</?, -<p), C : (z, </?) -> (2, <p). (51) 
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B Calculation of the renormalized scaling di- 
mensions 

In this appendix, we calculate the correction of scaling dimensions in Table 
1, up to the first order of y , in some degenerate cases. The derivation 



here is simpler than the original one |18], |J . Let us consider the following ID 
quantum Hamiltonian 

H = H + ^- [ L dvO u (52) 



2tt Jo 

where Ho is a fixed point Hamiltonian, L is the system size (in Table 1 I is 
related with L as / = log(L/a)), and Oi{= 0\) is a scaling operator whose 
scaling dimension is X\. We set the short-range cutoff as 1. According to 
Cardy [T9|] , the size dependence of excitation energies up to the first order 
perturbation is given by 

In 2n 
AE n (L) = — (x n + C n i n Xi(L) + ■•■) = T x n (L), (53) 

where x n is the scaling dimension of the operator O n , C n \ n is the operator 
product expansion (OPE) coefficient of operators O n and (D\ as 

1 (z, z)O n (0, 0) = C nln z- hl z- h *O n (p, 0) + ■ • ■ , 

in which h\(— h\) is the conformal weight of 0\ (x\ = 2h\). We used the 
notation 

Ai(l) = Ai (t) ' 

which comes from the renormalization group equation. 
For the sine-Gordon model, we denote 

x = ^( 1 + ^)- 

near K = 1 (m = 2) or K = 4 (m = 1). We can rewrite the Lagrangian 
density as 

C(z,z)=C (z,z) + C I (z,z), (54) 

where 

2 

Ttl 

C (z,z) = —(V<f ) ) 2 , 
16 



and 



and we set 



CAz, z) = —M + a/2 cos y/2m<> 

4vr 2y/2n 



Ao 



V2- 



(55) 
(56) 



B.l x p0 and x p3 

We calculate x p0 and x pl up to the first order of y's. We denote operators 
Vz cos m<\)j\pl (whose scaling dimension is x p0 ) and a/2 sin m0/ 'v2 (whose 
scaling dimension is x p3 ) as O c and (9 S respectively. We have the following 
OPE's 

M(z,z)O c>s (0,0) = ^O c>s + ---, (57) 
a/2 cos V2m<f){z, z)O c>s (0, 0) = ±^L— (0, 0) + • • • , (58) 

v 2 M 



where + is for O c and — is for O s . From these OPE's and eq. (p3|), we can 
obtain the scaling dimensions, 



V 7 2 2 2 a/2 a/2 V ; 



for a/2 cos m0/ a/2, and 



, P 3/r, _ 1 , 1 2/o(^) 1 y<p(L) 



x p-\L) = - + 1^11 - -^^^ (60) 
V ; 2 2 2 a/2 a/2 y J 

for a/2 sin m0/ a/2. Setting = yo(l + t), we have the scaling dimension 
described in Table 1. These are consistent with the results obtained by 



Giamarchi and Schulz |T8 



Similar calculation can be applied to x 



m5,m6,m7,m8 



B.2 Marginal operators 

For the BKT transition, there exists a hybridization between the marginal 
field Ai and the operator a/2cosa/277?.0 (m = 2 for K — 1, and m — 1 
for K = 4) ||. These two operators have the same scaling dimension and 
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symmetries at K = 1 or K = 4. Treating Lj as the perturbation term, we 
have the following OPE's 

C I (z,z)M(0,0) = ^-lV2cosV2m0(O,O) + ---, (61) 

2tt\z\ 2 

C I (z,z)y/2cosV2m(f>(0,0) = — ^v^cos V2m<p(0, 0) + ^r^M(0, 0) + • • • . 

2ir\z\ 27i\z\ 

Setting = y (l + 1), and diagonalizing these equations, we have the or- 
thogonal operators, 

- ~t) M + y| (l+ ~t) v^cos V2m0 (62) 
with the scaling dimension (up to the first order of y and t) 

x m0 (L) = 2 - y (L) (l + ~t 



3 



and 



^ (l - it) v^cos V2m0 - (l + ^) M (63) 
with the scaling dimension 

x ml = 2 + 2y (L) + 

These scaling dimensions are consistent with the results obtained by one of 
the authors |§. 

Scaling dimensions of remaining operators in Table 1 does not split. 

C Symmetry in twisted boundary conditions 

In eq. (|I7|), we showed the unitary transformation of TBC for spin operators. 
In this appendix we discuss the unitary transformation of other operators 
T R , P, T. They are not well defined except for the special $, in contrast to 
the spin operators. 
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C.l Translation operator 

The unitary operator U& of the twisted boundary condition (|17D is trans- 
formed with the translation as 

TrU^Tu]- 1 = exp^|:(j-l/2)^ +1 

= P$exp (-j; S t) exp^S 2 ), (64) 

therefore, we obtain 

U*T R [Us]- 1 = expiiQSm, (65) 

where we introduce the operator T R = exp(z$S'|i / L)Tr. The operator exp(z$5'j : 
is not well defined except for $ = 27rZ(Z:integer) 

U^TrIU^Y 1 = {-l) 2Sl T R . (66) 

In the notation T R , the periodicity of the energy and the momentum eigen- 
value under $ — > $ + 2it(S: integer) or $ — > $ + 4tt(S: half-odd-integer) 



44, 45 1 becomes apparent, in contrast to |46, 47]. Note that when we define 



the wavenumber q' = q + ^S^/L as 

T^^S^)) =exp{Kq + *Sl/L))\'t!>(q,SZ,$)), (67) 

Lq'/2ir is not integer in general, though it corresponds to the eigenvalue 
mn + n$/27r of the conformal spin operator L — L under TBC (see eqs. 
(|39|) , fl42"|) , dj3|) , (|50D ), and has a more natural symmetry structure than 



Tr, as will be seen later. 

C.2 Discrete symmetries P,T 

The unitary operator [/$ is transformed with the parity as 



PZ7*P = exp ya-l/fc. 



^ i=i 

= f7_*exp(i$S£), (68) 

therefore, we obtain 

C/$P[C/$] _1 = E/ 2 *Pexp(-i$S'f.), (69) 
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The operator exp(i<&Sif) is not well defined except for $ = 7rZ(/:integer) 

U^Pprt]- 1 = (-l) s T l U 2nl P, (70) 
where we use that Sj, is integer. In this case, we obtain 

(U 27Tl P) 2 = 1. (71) 



Similarly, we obtain 
and 



U^Tpt,}- 1 = U 2 *T, (72) 



(U 2 *T) 2 = 1. (73) 

For the twist angle $ = nl, Hamiltonian (|l^) commutes with the oper- 
ators U 2lT iP, U 2w iT. Moreover, these operators forms Z 2 group, that is, the 
eigenvalue of U 2w P, U 2n T is ±1. 

Next for the twist angle $ = nl, we discuss the relation between the 
operator U 2n iP and the translation operator T R . Using (0), we obtain 

T^V^PTl = exp(- r ^)T fl ^[T R ]- 1 T R PT /? 

= exp(2niS^l)U 2wl P (74) 

Since the operator exp(2iriSl) is ±7 for S integer or half-integer, this means 

T R (U 27Tl P) = (-l) 2Sl (U 27T P)(T R y\ (75) 

that is, in the momentum space for the even-integer 25/, the spectrum is 
symmetric with respect to q — 0, tc, whereas for the odd-integer 2SI, the 
spectrum is symmetric with respect to q — ±7r/2. For the half-odd integer 
S case, under TBC ($ = tc) the ground state (q = 0) is exactly degenerate 
with a state q = tt [BH], which is related to the Gaussian transition |26| . 



C.3 VBS states 

As an application of the previous subsections, we discuss the generalized Z 2 x 
Z 2 symmetries or valence bond solid (VBS) |27| under the twisted boundary 
condition. It was shown that under the twisted boundary conditions the 
quantum numbers P, T are the good quantum numbers characterizing the 
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generalized Z 2 x Z 2 symmetries [2B|. Here we generalize these results for the 
twisted boundary conditions with the translational invariant case. 

The spin variable can be represented by the Schwinger bosons as follows 

Sj = -{cLjdj — bj~bj), Sj~ = a~j"bj, = ajb^, (76) 

with the constraint that the boson occupation number at each site ajaj+b'jbj 
is 25. 

The VBS states with the TBC can be written as 
\S, M, TBC($ = tt)> = (a+bfe- i7T/2L - b + L a+e™ /2L ) s ~ M 

L/2-1 

x n K-i^ e " W2L - ^-i4^ /2i ) s+M 

3=1 

X ( a 2jHj+l e ln ^ L ~ ^2j a tj+l ei7T ^ 2L ) S 

x (ai_ 1 bte- i *' 2L -bi_ 1 ate**/ 2L ) s+M \Q) (77) 

where M is an integer for the integer S, or a half-odd integer fot the half-odd 
integer S (here we include bond-alternating cases). First we make a parity 
transformation for the VBS state 

P\3, M, TBC(<$> = tt)) = (-1) SL \S, M, TBC(<S> = -tt)) (78) 

where we use P\0) = |0). Then twisting with U^, we obtain 

U 27T P\S, M, TBC(<5> = tt)) = (-1) SL - S+M \S, M, TBC(<5> = tt)) (79) 

where we use ?72tt|0) = |0), and ^■n-aibtU^J' = a^bf exp(27ri(L — 1)/2L). The 
same discussion applies for \J 2 ^T. Therefore, each M-VBS states is charac- 
terized by the discrete quantum numbers U 2n P = U 2n T = (^—i^SL-s+M _ 

Similarly we can classify the intermediate large D phase with the discrete 
quantum number U 2n P, U 2n T under TBC. 
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Table 1: Operator content of the sine-Gordon model at K — 1 and K = 4. 
Here we consider the BKT transition of the y<f,{l) = Vo(l) branch, and we 
denote the deviation t from the BKT critical line as y<f>(l) = yo(0(l + (f° r 
the ^(/) = — £/o(0(l + branch, the role of operators x p0 <-> a; p3 , x m5 ' m6 <-> 
x m7,m8 interchanges). *:Strictly speaking, these operators (x m0 ,x ml ) are hy- 
bridized under renormalization. t : P,T should be interpreted as L^T, t/ 27r P 
under TBC. 

Operator in Operator in Renormalized Dicrete 

s-G model (K=l) s-G model (K=4) scaling dimension symmetries^ 

Oi, + 0-i,o O 1/2i0 + O_ 1/2>0 (TBC) ± + |y o (0(l + §*) P = T = 1 

* pl ' p2 O 0)± i Oo,±2 3-30»(O P = 1 

* p3 Oi,o- O-i,0 O 1/2;0 - 0_ lA o(TBC) | - isrt>(0(l + 2*) P = T = -1 



x cl ' c2 ±1 , ±1 0±i/2,±2(TBC) 1 

x m0 * ~^0(t)d(j) -*d<i>d<i> 2-y (0( 1 + |*) P = T = l 

x mi* O 2 ,o + O_ 2 ,o Oi >0 + O_i )0 2 + 2y (/)(l + |i) P = T=1 

x m2 O 2 ,o-O_ 2 , Oi,o- O-i,0 2 + y (0 P = T = —1 

x m3,m4 2 - y (l) P=l 

x m5 ' m6 ao±i,±i + aOzpi.ii ao ±1A±2 + ao TlA±2 (TBC) 2 + yo (0(i + *) p = 1 

x m7, m8 oo ±l±l _ ao T i,±i ao±i /2 ,± 2 - ao T i A j 2 (TBC) 2- yo (/)(i + t) p = -i 



excitations 


dA£/dA| A=0 


expected ratio 


pl,p3 


0.50426315 


1/2 


m3,m0 


1.23831439 


4/3 


m7,m0 


0.30949032 


1/3 


m3,m7 


0.92882406 


1 


m5,m2 


0.93899885 


1 



Table 2: dA£/dA| A =o for L = 16 system. 
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Figure 1: Excitation energies of = ±2, q = 0, P = 1 state with PBC 
(x), and of S t z ot = 0,q = states with TBC (• : U 2w P = U 2w T = -1 and 
o:U 2w P = U 27T T = l). 
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Figure 2: Excitation energies of = 0, q = 0, P = T = 1 states with PBC 
(o), of S% = 0, q = 0, P = T = -1 state with PBC (•), of = ±2, q = 
states with TBC (x), and of S% = ±4, q = 0, P = 1 state with PBC (+), 
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Figure 3: Energy difference E c0 — E cl for L = 16 systems. 
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Figure 4: Size dependence of the averaged scaling dimension (a: p0 -|-3x pl ' p2 )/3. 
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Figure 5: Excitation energies of L = 12, D = near the XY-Haldane tran- 
sition point, x's are Sj> = ±2, q — 0, P — 1 excitations under PBC, 
o's are = 0, U^P = U2 n T = 1 excitations under TBC, and »'s are 
S% = 0, U 2 nP = U 2 nT = -1 under TBC. 
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Figure 6: Size dependence of the crossing point. The extrapolated value is 
A c = 0.966, which is the transition point between the S = 2 Haldane gap 
and the XY phases. 
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Figure 7: Size dependence of the averaged scaling dimension (x p0 + 3a; pl )/4 
for A r 
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Figure 8: Excitation energies of L — 10, A = 1 near the Haldane-XY and 
XY-large D transition points, x's are = ±2, q — 0, P — 1 excitations 
under PBC, o's are = 0, U 27r P = U 27T T = 1 excitations under TBC, and 
• 's are = 0, U 27T P = U 27r T = -1 under TBC. 
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Figure 10: Size dependence of the crossing point. The extrapolated value is 
D C 2 = 2.39, which is the transition point between the XY and the large D 
phases. 31 
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Figure 11: Size dependence of the averaged scaling dimension for D cl . 
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Figure 12: Size dependence of the averaged scaling dimension for D c2 
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Figure 13: Size dependence of the Gaussian fixed points in XY phase on the 
A = 1 line 
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